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Abstract. We show that there exists a choice of gauge in which the electromagnetic 4-potential may be written as the 
difference of two 4-velocity vector fields describing the motion of a two-component space-filling relativistic fluid. Maxwell's 
equations are satisfied immediately, while the Lorentz force equation follows from the interactions of sources and sinks. The 
usual electromagnetic quantities then admit new interpretations as functions of the local 4-velocities. Electromagnetic waves 
are found to be described by oscillations of the underlying medium which can therefore be identified with the 'luminiferous 
aether'. The formulation of electrodynamics in terms of 4-velocities is more general than that of the standard 4-potential in 
that it also allows for a classical description of a large class of vacuum energy configurations. Treated as a self-gravitating 
fluid, the model can be explicitly identified with Nelson's stochastic formulation of quantum mechanics, making it a promising 
candidate as the classical field theory unifying gravitation, electromagnetism and quantum theory which Einstein had sought. 



INTRODUCTION 

Motivated by the fact that macroscopic waves tend to propagate in a medium of some kind, many physicists have in 
the past attempted to find a description of electrodynamics in which electromagnetic waves may also be described 
by the motion of an underlying medium, commonly referred to as the 'aether' Oil. Despite significant efforts, no such 
description was found, and the negative results of the Michelson-Morley experiment|2] put the final nail in the coffin 
of the aether concept as it was then understood. The failure of this program eventually led to the introduction of the 
concept of a 'field' |31 requiring no underlying medium, and on quantisation of these fields, to 'quantum field theory'. 

We demonstrate here that Maxwell's equations can in fact be derived from the motion of an underlying medium if we 
assume from the outset that the underlying spacetime is Lorentzian rather than Galilean. While our formulation differs 
from earlier conceptions of the aether, it does show that an alternative description of electrodynamics did (and does) 
exist which did not require the introduction of the field concept. The medium can be interpreted as a two-component 
relativistic fluid with sources and sinks playing the role of charges. This fluid dynamical model allows an explicit 
connection to be made with Nelson's stochastic formulation of quantum mechanics, suggesting that the foundations 
of quantum electrodynamics may lie purely in classical gravity as Einstein had long believed. 

We assume a metric with signature (+, — , — ,—), and follow the conventions of Jackson[4J throughout. 



THE RELATIVISTIC CONTINUUM 

We first establish the reference system and coordinates that will be used to describe the motion of the continuum. We 
then define the 'continuum gauge', showing how Maxwell's equations equations appear. We also show that the Lorentz 
force equation translates into a third order partial differential equation constraining the continuum's motion. 

Coordinates and Reference Frames 

Consider an arbitrary relativistic inertial frame with 4-coordinates = (cf ,jc,y,z), so that the spacetime partial 
derivatives are given by = (7^,— V) and d^^ = (^r^jV) respectively. Suppose that there exists a continuous 



medium in relative motion to this frame which spans the entire spacetime. To describe positive and negative charge 
configurations we will see that this continuum must consist of two superimposed but mutually independent components 
which we shall refer to as the 'positive continuum' and the 'negative continuum' respectively. 

Let T be the proper time in the inertial frame, and consider the instantaneous motion at proper time T of either 
component of the continuum at any 3-position r — (x,y,z)- Then the 3-velocity of that component at r as measured by 
the inertial frame is, 

where t is the time as measured by a clock moving with the continuum component. We can therefore define the interval, 

ds^ = c-dx^ = c^dt^ - dx^ - dy^ - dz^ . (2) 
Similarly, we can define a 4- velocity vector field describing the motion of the continuum component as, 

,, dx^^ f dt dr\ , . 

where 7= (1 — v^/c^)^^/^ is the Lorentz factor at each point. This 4-velocity clearly satisfies, 

u^'u^=c^, uf^'^u^^Q, (4) 

where partial derivatives d^u^ are written as u^'^ for convenience. 

The above definitions and identities hold for both the positive and negative continua independently, and henceforth 
we distinguish the two sets of variables by a '+' or '-' subscript respectively as necessary. 

The Continuum Gauge 

The key to what follows is to split the electromagnetic potential 4-vector A'' into the sum of two components 
and , which we identify (up to a dimensional constant k) with the two continuum 4- velocities u'^ and u^, 

A^=A'i^+A^, where A^ ^ ku^ ^ {(!)+, A+) , A^ ^ -ku^ ^ {^^,A^) . (5) 
The condition (0) implies the following covariant constraint for both A'I^ and A^l, 

A'i^A+f,=A'iA-^=k^c^ . (6) 

We will refer to conditions (|5j and ^ as the 'continuum gauge'. This is a non-standard choice of gauge, and we 
will demonstrate its consistency in the next section. The antisymmetric field-strength tensor can now be defined as, 

P^iv ^^v.ti_^ti.v _(E,B). (7) 

Other standard properties now follow in the usual way. From the definition Q, F^^^ satisfies the Jacobi identity, 

ptiV,l_^pvX.^i^pXti.v (8) 

and this is just the covariant form of the homogeneous Maxwell's equations. One can define the 4-current as the 
4-divergence of the field-strength tensor, 

J^'-i^Fy =(^PJ), (9) 

and this is the covariant form of the inhomogeneous Maxwell's equations. Charge conservation is guaranteed by the 
antisymmetry oiF^^ . The covariant Lorentz force equation takes the following form, 

^ = (10) 



where Q, M and V — (cjv, }VV) are the charge, mass and 4-velocity vector of the observed particles. This cannot be 
derived directly from the definition of the 4-potential, and must be considered for now as an auxiliary constraint. 
The charge 4-velocity V and scalar charge Q are related to the 4-current density through the following equation, 

J^ = QV^, where y^'y^:=c^ V°>c. (11) 
The constraint on allows us to separate the 4-current uniquely into the charge and its 4-velocity. Indeed we have, 

/ 1 \ 1/2 

e = sgn(/) • f -^7% j , (12) 

where the sign of the 0-component of the 4-current appears to ensure that the 0-component V" of the charge 4- 
velocity is positive. Since the sign of 7" cannot be flipped by a Lorentz transformation, each 4-velocity vector field 
can only account for either positive or negative charge configurations, and hence the need for two separate continuum 
components, each with opposite contributions to the 4-potential. The gauge based upon a single 4-vector field was 
precisely that introduced by Dirac in his classical model of the electron! 5], and it is noteworthy that he was also led 
to speculate that this 4-velocity field described the motion of a real, physical, 'aether' [6]. The form of the charge 
4-velocity in terms of the continuum 4-velocity now follows directly from il l\ . 

Besides the mass M which is determined by initial conditions, each of the terms in ilO\ may be written in terms of 
the 4-velocities u'^ and u^. From the definitions of F^^, J^, Q and V^, we find that the Lorentz force equation dlOt 
translates into a complicated third order partial differential equation constraining the 4-velocities. The conservation of 
mass follows from the continuity equation for mass density, 

{MVn.n^O, (13) 

which is ensured if the flow of mass density follows the flow of charge density. We will see later that the Lorentz force 
equation follows from the fluid dynamical interactions between sources and sinks, and this will complete our picture 
of classical electrodynamics in this gauge. 



THE CONSISTENCY OF THE CONTINUUM GAUGE 

We have identified the components and A'^ of the 4-potential with the 4-velocities u'l and u'^ of the continuum 
satisfying the conditions and and have referred to this gauge choice as the 'continuum gauge'. It is not obvious 
that this gauge choice can be applied consistently to all electromagnetic field configurations, so we demonstrate its 
consistency here, and give explicit solutions given for the point charge and the plane electromagnetic wave. 



Proof of Consistency 

In order to prove consistency, it is necessary to find a decomposition of the 4-potential as the difference of two 
4-velocity fields satisfying equations Q and simultaneously. Using the notation of Q, we therefore need to find, 
given any 4-potential A'' = {(p,A) defined up to a gauge transformation A'' ^ A'' +(9''v/, two 3-velocity fields v+ and 
v_ satisfying the following conditions, 

A 

^ = 7+ -7-, = 7+V+ - 7-V- . (14) 

The second of these equations is a simple geometrical vector identity, and it is clear that any solution set for 
(7+v+,7-V_) will form a surface of revolution about the axis defined by A. To find the solution surface explicitly 
for a given (0, A), it is convenient to take the origin to lie at A/2k, and to use polar coordinates (r, 9) in any plane 
containing A, where r £ [0,°°] is the radial distance from the origin and 6 e [0, n] is the angle made with respect to the 
direction of A. Note the following simple chain of identities. 




:yf^^7=\/l + (^^)', (15) 



so that from ([T4l we have, 

Applying standard trigonometric identities to our geometrical picture, we obtain, 

(r+v+f = r^+A^ /4-k^ + ^ cose, (7_v_)2 =r2+AV4/t2-ycos0, (17) 
so that the set of solutions on the plane in question is determined by the condition. 



(j) ^ ^A^/4+Akrcose+k^-{r^-+c^)-^A^/4-Akrcos9+k^{r^+c^). (18) 

Note that given any solution for {(p,A), a solution for (— ^,A) is obtained by letting 9 ^ K — 9. Note also (i) that 
(p =0 whenever 9 = n/2 including when r = 0, (ii) that for a given value of r the magnitude of (p is maximum when 
= 0, (iii) that for 9 — 0, (p is a monotonically increasing function of r, and (iv) that (p ^ Acos 9 as r ^ °o. 

In conclusion, for a given value of A = | A|, equations il4l will have solutions whenever 1 1 < A. In the special case 
^ = the solution surface for 7+v+ is just the plane perpendicular to A passing through the point A/2k, throughout 
which |v+| = |v_|, and |7+v+| >A/2k. For other values of \(p\ < A the solutions form a paraboloid-like surface of 
revolution about the A axis. The sign of (j) determines which side of the 9 = ±7z/2 plane the solution surface lies. 

It is always possible to choose the function Xj/ defining the choice of gauge in such a way that (f) —0 every whereQ. 
Since solutions to J14> always exist in this case, this proves that the continuum gauge is indeed a consistent one. 

It is important to note that there is actually a significant additional degree of freedom inherent in the way the 
decomposition of A^ is made into 4-velocity fields, which goes beyond the standard gauge freedom. First of all, for 
each electromagnetic configuration there will be a continuum of gauge choices for which a continuum gauge solution 
set exists. Secondly, for any particular choice of gauge for which a solution does exist, there will in general be an 
entire two-parameter surface of possible solutions for v+ and v_ at each point in space. We will show later that these 
velocity vector fields correspond to the motion of massive discrete particles, so that this freedom may have a real 
physical significance as a possible classical source of dark matter. 



The Point Charge 

Let us now find the vacuum configuration which describes a positive charge q positioned at the origin. The 
corresponding electromagnetic fields are given by, 

E=^ = -V(i), B = 0. (19) 
We seek a 4-potential of the following form which only has contributions from the motion of the positive continuum, 

A^+ = i^+,A+) = {kcy^kyy) , A^_ = ((/)_,A_) = {-kc,0) , (20) 
where the velocity vector field v is to be found. The corresponding electromagnetic fields E and B are given by, 

E = = -V{kcY)--^{kYy), B = VxA+ = Vx(;t7v). (21) 

c at c at 

For any electrostatic configuration with stationary charges we have B = V x [kyy) = 0, so there must exist a scalar 
field i/A such that ky\ = Vy/. After some algebraic manipulation this can be seen to imply that, 

so that in terms of y/, the E field is given by, 

E = - V ({kV + (Vv/)2) '/^) (Vv/) . (23) 



c 



Because of the rotational and time invariance of the problem, we need only look for solutions of the form \ir = y/{r), 
so that Vxj/ = dij/fdr and the second term of (I23> vanishes. Comparing with ( I19t . it is clear that Xj/ must satisfy, 




- + «, (24) 
r 

where a is a constant of integration. Since the charge is positive and the velocity of the continuum should vanish at 
infinity, we require a = kc for a real solution to exist. From the resulting differential equation for y/ is as follows. 



where either the positive or negative square root may be chosen, as the 4-potential depends only on the magnitude of 
the velocity and not its direction. There is therefore insufficient information to specify whether the positive charge acts 
as a source or a sink (or both). The solution for the velocity field and the corresponding Lorentz factor is therefore. 



V 



1/2 



- = ±(i-(i + f)")"'", r-i + f. (26) 

c \ \ krc/ J krc 

Note that q/krc becomes singular at the origin, implying that the continuum velocity in (I26> becomes equal to c there. 

The above confirms that the electromagnetic fields outside a positive point charge can indeed be described by the 
motion of the positive continuum, and that the corresponding potential 4-vector A'^ is expressible in terms of the 
4-velocity u^. An identical calculation can be performed to show that an analogous result is true for negative charges. 



The Plane Electromagnetic Wave 

While in principle one can claim that all electromagnetic configurations ultimately originate from the presence of 
charges, there do exist nontrivial configurations in which no charges are present, the most obvious and important 
example being that of the electromagnetic wave. It is therefore important, both for this reason and from a historical 
perspective, to show explicitly how plane waves arise in the present context from the motion of the relativistic 
continuum. We turn to this problem now. 

Let us consider a plane electromagnetic wave with wave-vector k travelling in the x-direction with the E-field 
plane-polarised in the y-direction. The 4-potential describing this plane wave is, 

= (0,A) = (0,0,A,.cos(a)f- K-x),0), (27) 
(where CO —ck), with corresponding E and B fields, 

E = (0,£v,0) = (0,K-AvSin(a)r-K-x),0), B = (0,0,5^) = {0,0,KAySin{m - kx)) . (28) 

We therefore seek solutions of the form, 

Al = {kcY+,kr+y+) , A>i = {-kcY--kY^y-) . (29) 

Applying and equating with j27l l we obtain the two conditions, 

Y+ = Y-, ^(7+v+ - 7-V_) = {Q,AyCos{cot - kx),Q) . (30) 

Ignoring equal velocity motions of the positive and negative continua which have already been shown to have no 
electromagnetic consequences, these conditions allow us to restrict our attention to solutions of the form, 

v+ = -V- = (0,v,0), where -= , / ^ . (31) 

c VA2 + 4c2 

and we have defined A = A,, cos (of — Kx) for convenience. The velocities of the positive continuum and the negative 
continuum here are equal in magnitude and opposite in direction, so that there is no net charge, with the motion of both 
being parallel to the electric field but =F^/2 radians out of phase respectively. It also follows from ( I31> that the velocity 



of the continuum can never exceed the speed of Hght, irrespective of the intensity of the plane wave. Substituting ( 13 U 
into i29\ the motion of the continuum is given by, 

u1:^{^c^/k^+A2/4k^,0,A/2k,0), = {^c^/k^+A^/4k^,0,~A/2k,Q). (32) 

These equations clearly show that the propagation of a plane electromagnetic wave is described by the oscillation 
of the medium in the direction of the electric field - the positive continuum oscillates 7z/2 out of phase with E while 
the negative continuum oscillates with the same magnitude and precisely the opposite phase. Thus the propagation of 
electromagnetic waves is seen to be a direct manifestation of the oscillations of the underlying relativistic continuum. 



Gauge Redundancies and the Principle of Superposition 

While the usual principle of superposition obviously still holds for the 4-potential, we can now supplement this with 
the following continuum-gauge-inspired superposition principle. 

Consider two 4-potential fields A^^ = (c7+ — c7_,7+v+ — 7-V_) and A'^^ — (c7^ — cy_,'y'_^\'_^ — 7^v'_) in the 
continuum gauge which describe two different 4-velocity field configurations. Then the superposition of the two field 
configurations is described by the 4-potential A"^ = (c7^ — cy^,y_l\'l_ — yW'_) where the velocity vector field v" 
(respectively v" ) is given by the pointwise relativistic sum of v+ and v'^ (respectively v_ and v'_), 

" (33) 



± l+V±-v'±/c2- 

As mentioned earlier, the description of an electromagnetic configuration in terms of 4-velocities and is far 
from unique, as for each of the infinite number of 4-potentials A** ~ (0,A) with < |A| describing that particular 
configuration, there exists an entire two-parameter set of solutions at each point. 

Recall the particular gauge choice in which 0=0 everywhere. We saw that the simplest 'lowest energy' solution is 
given in this case by v+ — — v_ = A/2k. However, we also saw that it is possible to add, relativistically in the sense 
of j33> , the same, arbitrary, possibly time-dependent, 3-velocity vector field to both v+ and v_ without changing the 
4-potential. If these velocity fields have a real physical meaning then this additional freedom will correspond to a large 
class of vacuum configurations which can perhaps be interpreted in terms of the motion of an arbitrarily distributed 
'Dirac sea' of particles and antiparticles. This provides a means of adding energy density to the vacuum without any 
observable electromagnetic effects. 



THE CONTINUUM AS A RELATIVISTIC FLUID 

In this section we show that the spacetime continuum must be a relativistic fluid of massive discrete particles, and that 
interactions between sources and sinks give rise to the Lorentz force equation. The fact that both Maxwell's equations 
and the Lorentz force are consequences of the relativistic fluid model is a strong indication that there is more to this 
description than mere formalism, and that classical electrodynamics may in reality have a fluid dynamical basis. 



The Massive Continuum 

We saw in i26i that the velocity of the continuum decreases with radius outside of the point charge acting as its 
source. Had the continuum been massless, its velocity would have been constant and equal to c everywhere. We 
therefore conclude that the continuum has mass and that there is an attractive central force acting on the continuum 
outside of the charge. 

It is possible to derive an expression for this attractive central force. In particular, if we assume the charged particle 
is centred at the origin, then the force f acting on an infinitesimal element of the continuum at radius r must satisfy [7J, 



(34) 



where m — p,nSV is the mass of the test element assuming that it has mass density p,„ and occupies volume 5V. To 
find the value of dv/dt, solve J26l l for and differentiate the resulting equation with respect to t to find an expression 
for dv/dt in terms of v. Rearranging terms and simplifying, the field at radius r is found to have the form. 

Thus there appears to be a Coulombic attraction between the charge and the continuum around it, with the continuum 
having a charge-to-mass ratio of —c/k. This is quite mysterious as in our model charge is defined in terms of the motion 
of the continuum, so clearly the continuum itself cannot be charged. The mystery will be resolved in due course. 

Assuming continuum conservation, the continuum density p„ will satisfy the following continuity equation. 



5m(P«^'') = %^-V.(p„7v)=0 ^ ^|-(r2p„7v)=0. (36) 
at or 



where we ignore the time-derivative term as the system is in a steady state condition, and use the rotational symmetry 
to rewrite the divergence term in its spherical polar form. The solution is, 

where 5 is a radius-independent proportionality factor. Now, the flux of continuum passing through a spherical shell 
at radius r is just <I> — Anr^PnYv (where the factor of 7 takes into account to the relativistic contraction in the radial 
direction). But this is precisely the constant S in (I37> which can therefore be identified as the strength of the charged 
particle sink/source. 



The Relativistic Fluid 

We discovered in the previous subsection that there is an inverse-square law attraction of elements of the continuum 
towards the point charge. Given that the continuum density is greater closer to the charged particle, let us investigate 
the possibility that the continuum may be a continuous, compressible, medium whose attractive self-interactions result 
in the observed attraction. If the attractive force between two volume elements of the continuum is given by, 

dF {ri,r2)^ Pn (r 1 )p„ (r2)/( |r 1 - r2 1 ) . (38) 

where f{r) is some polynomial in r, then a little calculation shows that an inverse square attraction is possible only if 
f{r) ^ r^'^. However, the magnitude of the resultant force on any element turns out to be infinitely large in this case. 

There are three sources of these (logarithmic) divergences - (i) the contribution from the core of the point charge 
where the continuum density becomes infinite, (ii) the contribution from the continuum at infinity, and (iii) the 
contribution from continuum elements in the immediate neighbourhood of that element. The first of these can be 
avoided if the charges are not pointlike, the second can be avoided if the universe is either bounded or homogeneous, 
and the third can be avoided by discarding the idea that the continuum is some kind of continuous elastic medium, but 
rather consists of a fluid of interacting discrete particles. 

We are therefore led to conclude that our relativistic continuum is a space-filling relativistic fluid and the electro- 
magnetic 4-potential must be defined in terms of the ensemble motion of the fluid as opposed to the motion of the 
individual discrete particles. If the instantaneous fluid velocity ar is (x), then the 4-velocity appearing in (|5} is, 

u^'{x)=<C^>, (39) 

where < > indicates the time-averaged motion of the particles in the neighbourhood of . All other electrodynamic 
quantities must be defined as time-averages in the same way. 

Although the configuration representing a charged particle is in steady state, the fluid itself remains in constant 
motion. Recall that the motion of an individual particle in the co-moving frame is described by the total derivative|8|, 

^ = (C'^v)^ - iCvd'K" - l:d^iCCv) = -{d^C - d'CnCv , (40) 
dz 2 

where we have added a vanishing term using the fact fliatC^C =c^-If 

we now consider the time-averaged version of 

( l40t and recaU flie definitions (|3|}, (|5} and (|7}, we find that. 



k—^^F^^Uy, (41) 

dx 

which is in the form of the Lorentz force equation. In particular we find that, on average, each particle moves as if 
it were charged with qjm — —c/k. This is precisely the charge-to-mass ratio observed in the Coulomb-like attraction 
of ( I35> . and so the earlier mystery has been resolved. Because ( I40> is a basic identity valid for any motion of the 
relativistic fluid, this conclusion holds irrespective of the precise nature of the interactions between the fluid particles. 

As further compelling evidence that classical electrodynamics has relativistic fluid dynamics as its basis, we will 
now show that the Lorentz force equation emerges automatically from the interaction between sources and sinks when 
they are not assumed to be fixed in position. 

The integral momentum equation for a fluid tells us that the force on a target charged particle with charge Q' due 
to a source particle of charge Q at distance r is given by the rate of change of momentum transfer to the target by the 
particles entering or leaving the source. If we suppose that the target particle has an effective radius R then, assuming 
spherical symmetry, it will have an effective volume of ^kE?. In accordance with (I37> . the density of fluid particles 
encountering the tai-get at distance r from the source is p„(r). If we further assume that each fluid particle is identical 
with mass m, then the 3-momentum carried by each is given by myv. Finally, the collision rate will be determined by 
the strength S' of the target. Thus the force on the target will be given by the product of these contributions. 



F - ^ KR^p„myvS' = , (42) 

where he have used ( I37> . This takes precisely the form of Coulomb's law if we make the following identification. 



where the charge Q is expressed in terms of the strength of the source S, the mass m of the fluid particles and the 
effective charge radius R. Clearly for \A2\ to hold, positive charges must effectively act as sinks, and negative charges 
as sources, or vice versa'. The validity of Coulomb's law in turn impUes the validity of the Lorentz force equation! 9], 
as we have assumed from the outset that relativity holds. This completes our description of classical electrodynamics. 



Stochastic Quantum Mechanics from the Self-Gravitating Fluid 

We have shown that classical electrodynamics can be explained at the macroscopic level in terms of the ensemble 
motion of a relativistic fluid of massive (i.e. gravitationally interacting), discrete particles. At a microscopic level this 
allows us to make an explicit identification of our model with the stochastic formulation of quantum mechanics. 

In his monograph! 10!, Nelson gave a detailed derivation of quantum mechanics on the basis of the conservative 
diffusion of a classical fluid, wherein the Schrodinger wavefunction is identified with the density of the fluid thus, 

V/=Vpe'^/^ (44) 
where S is the stochastic analogue of Hamilton's principle function. However there remained a number of important 
unresolved problems. The first of these was "...to find a classical Lagrangian, of system + background field oscillators 
+ interaction, that [...] produces a conservative diffusion system. " Traditionally the system is assumed to be coupled 
in some way to an electromagnetic background! 1 1 !• We have succeeded here in showing that the classical fluid is the 
electromagnetic background which Nelson sought. 

Neither was he able to come to any definite conclusion about the nature of the interparticle interactions responsible 
for the conservative diffusion - except that it could not be gravitational on dimensional grounds and that it was possibly 
of electromagnetic origin. This second issue has also been resolved here as gravitational interactions are themselves 
responsible for electromagnetism. The possibility of a gravitational explanation for the conservative diffusion, and 
consequently for the observed magnitude of Planck's constant, has also been proposed by Calogero ltl2il . 



The fact that equal velocity contributions of fluid particles from the positive and negative continua have no electromagnetic effects means that the 
net momentum transfer must be zero, which in turn implies that particles in the negative continuum must have equal and opposite mass to those in 
the positive continuum. The constant in JsJ must then be proportional to the mass of the fluid particles, so that the 4-potential A'' is nothing but 
the net 4-momentum of the two-component fluid. This is a radically different interpretation of the 4-potential from the one we are accustomed to. 



SUMMARY AND CONCLUSION 



We have demonstrated the simple yet profound result that all of the equations of classical electrodynamics follow 
from the motion of a two-component relativistic continuum satisfying the standard equations of relativistic fluid 
dynamics. Charged particles appear as sources and sinks of the continuum in this framework, while electromagnetic 
waves are associated with oscillations of the continuum. There is a freedom inherent in the 4-velocity description 
of electrodynamics which can potentially account for 'dark matter'. The identification of the vacuum as a self- 
gravitating fluid of discrete particles makes possible an explicit connection to Nelson's stochastic formulation of 
quantum mechanics, raising the tantalising prospect of having a natural, unified, description of gravitation and quantum 
electrodynamics purely in terms of classical general relativity. 

Maxwell and others had struggled to find a mathematical description of the underlying medium, the 'aether', in 
which electromagnetic waves were presumed to propagate. Although the continuum we have described is not precisely 
equivalent to the notion which the earlier proponents had had in mind, our analysis does show that a description of 
electrodynamics in terms of an underlying continuum is possible. This is particularly important as the failure to find 
such a formulation historically contributed to the origin of the concept of 'fields' postulated not to require such a 
medium. The field concept may not have been necessary after all. 

A number of issues still remain open. We are still left to ponder the existence, interpretation and physical properties 
of the fluid particles and their sources and sinks, to explain the origin of quantised mass and charge, and the presence 
of the two continuum components. Considerable evidence has already been gathered that each of these issues can 
be resolved completely within the framework of general relativity, and we feel that a unified classical description of 
quantum theory and gravity is now close at hand. It would certainly be fitting if Einstein's dream were finally to be 
realised on the 100th anniversary of the birth of his theory of relativity. 
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